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Abstract 

A formalism for target fragment transport is pre- 
sented with application to energy-loss spectra in thin 
silicon devices. Predicted results are compared with 
experiments using the surface-barrier detectors de- 
veloped by J. P. McNulty. The intranuclear-cascade, 
nuclear- react ion model does not predict the McNulty 
experimental data for the highest energy events. A 
semiempirical nuclear cross section gives an ade- 
quate explanation of McNulty’s experiments. Appli- 
cation of the formalism to specific electronic devices 
is discussed. 

Introduction 

The early suggestion that some spacecraft anom- 
alies may result from the passage of the galactic ions 
through microelectronic circuits (ref. 1) has now been 
well established. Although the direct ionization by 
protons appears as an unlikely candidate, the recoil 
energy of nuclear-reaction products is a suspected 
as a source of single-event upset (SEU) phenomena 
(refs. 2 to 4). As a result, a number of fundamen- 
tal experimental and theoretical studies were under- 
taken to better understand the phenomena. McNulty 
and coworkers examined the energy deposition of pro- 
ton reaction products in Si using surface-barrier de- 
tectors of various thicknesses from 2.5 to 200 fi m 
(ref. 5). They also developed a Monte Carlo code 
for theoretical evaluation of energy deposition from 
such products. (See refs. 5 and 6.) A comparison of 
McNulty’s work with the well-established medium 
energy cascade code (MECC-7) developed by Bertini 
and coworkers at Oak Ridge National Laboratory 
showed some differences in predicted reaction prod- 
ucts and even greater differences in energy spectral 
contribution. (See ref. 7.) An evaluation of Si re- 
action products was likewise made by Petersen 
(ref. 4), and, although no direct comparison was 
made with McNulty’s experiments, an estimate of 
SEU rates in the trapped-proton environment was 
made. 

Following these fundamental studies were more 
detailed applications to specific-device geometries 
and parameters. Bradford evolved an energy depo- 
sition formalism (ref. 8) using the cross sections of 
Hamm et al. (ref. 7). McNulty et al. (ref. 5) ap- 
plied their Monte Carlo model to dynamic random 
access memory (DRAM) devices with reasonable suc- 
cess and discussed the implications of heavy- ion SEU 
phenomena on proton-induced SEU events through 
secondary reaction processes (ref. 9). The funda- 
mental consideration is the evaluation of the energy 
deposited within the sensitive volume (depletion re- 
gion) of the device in question as the result of a pass- 


ing proton. The ionization due to the proton itself 
makes only a small contribution to the critical charge. 
Nuclear-reaction events usually produce several re- 
action products (a heavy fragment and several lighter 
particles, although a few heavy fragments may be 
produced simultaneously on some occasions), and all 
the resultant products can make important contri- 
butions to the deposited energy. Such nuclear-event 
products are, of course, correlated in both time and 
space. 

There are three distinct approaches to a funda- 
mental description of the energy deposition events. 
McNulty and coworkers developed a Monte Carlo 
code in which multiparticle events are calculated ex- 
plicitly, including spatial and specific-event (tempo- 
ral) correlation effects. Although this is the most 
straightforward way of treating the full detail, it is 
a complex computational task. A second class of 
methods begins with the volumetric source of colli- 
sion events and calculates the SEU probability using 
the chord-length distribution. (See refs. 8 and 10.) 
Although in principle the correlation effects could be 
so incorporated, they appear to be ignored in both 
the cited references. A third approach in which linear 
energy transfer (LET) distributions and chord-length 
distributions are used seems most appropriate for ex- 
ternal sources. (See refs. 11 and 12.) This last ap- 
proach applies if the LET distribution from external 
sources is constant over the sensitive volume, but its 
applicability to volumetric sources is questionable. 
At the very least, this approach ignores correlation 
effects. 

Nuclear data bases for biological systems were 
examined in reference 13. The MECC-7 results 
underestimated by nearly a factor of 2 the energy- 
transfer cross section for multiple- charged ion prod- 
ucts. In a more detailed analysis (ref. 14), the 
Silberberg-Tsao fragmentation parameters were found 
to be superior to the MECC-7 results. The primary 
differences appear for the lighter of the multiple- 
charged fragments. Further comparison with experi- 
ments on Al targets shows both Monte Carlo nuclear 
models (McNulty’s code OMNI as well as MECC-7) 
to underestimate production cross sections for prod- 
ucts lighter than fluorine in proton- induced reactions. 
Although these intranuclear-cascade models are ca- 
pable of representing multiparticle correlation, the 
inherent inaccuracies in predicting cross sections is a 
serious limitation. 

In the present paper, the effects of nuclear recoil 
on electronic devices are examined and the develop- 
ment of a formalism for application to specific-device 
parameters is begun. As a test of our methods as 
they develop, the results are compared with the ex- 
perimental measurements of McNulty et al. (ref. 5). 



Symbols 


A 

atomic mass, amu 

Ap 

fragment atomic mass, amu 

CMOS 

complementary metal-oxide 
semiconductor 

dF 

de 

total absorption spectrum, MeV -1 

da 

dA 

cross section, mb 

E 0 

average fragment energy, MeV 

F z 

inward-directed flux at boundary of 
type-z ions 

f(Eo) 

spectrum average energy, MeV 

I 2 L 

integrate- integrate logic 

L 

feature size, fim 

MECC-7 

medium energy cascade code 

NMOS 

N-channel metal-oxide semiconductor 

OMNI 

McNulty code 

PMOS 

P-channel metal-oxide semiconductor 

P{E) 

energy-loss spectrum, MeV 

Qc 

critical charge, pC 

SOS 

silicon on sapphire 

z 

charge number 

€ 

energy loss, MeV 

a F 

nuclear fragment cross section, mb 


Microelectronic Upsets 

An electronic device is sensitive to the sudden 
introduction of charge into the active elements of its 
circuits. The amount of such charge that is sufficient 
to change state in a logic circuit is called the critical 
charge. As shown in figure 1, there is a rough 
relationship between the critical charge Q c and the 
device feature size L (ref. 11). This relationship is as 
follows: 

Q c « 0.0156L 2 (1) 

where Q c is measured in pC and L is measured in 
/im. Upsets in a device are then dependent on the 
charge produced in comparison to the critical charge. 

The charge released AQ in a material because of 
the passage of an energetic ion is related to the kinetic 
energy lost A E during the passage and is given by 


A Q = 


A E 
2243 


(2) 


where A Q has units pC and A E has units MeV. The 
energy lost by an ion in passing through a region is 

related to its stopping power = S z (E')j in the 

medium. The distance traveled before coming to rest 
is 

rE dpt 

If an ion is known to come to rest in distance x, then 
its energy is found through the inverse of relation (3) 
as 

E = R~ z l {x) (4) 

Equation (4) is used to calculate energy loss. The 
energy loss by an ion of charge Z and energy E in 
passing through the active region of a device with 
collection length L c is given by 

AE = E - R^[R Z (E) - L c ] (5) 

where 

L c = W e pi + W„ (6) 

In equation (6), is the epitaxial layer thickness 
and W n is the width of the depletion region (ref. 15). 
The energy loss depends on the particle isotope (i.e., 
ion mass) and angle of incidence. The ranga-energy 
relations described in reference 14 are utilized. As a 
practical matter to reduce numerical error inherent 
to numerical interpolation, 

A E = R Z 1 \R Z (E)]-R Z 1 [R Z (E)-L C ] (7) 

is used in place of equation (5). The result of 
equation (5) depends on the global error (fixed at 
1 percent) in the computer code, while equation (7) 
depends only on the local relative error (quite small). 
The charge introduced into the feature is given by 
equations (2) and (7). An example for a particular 
collection length of 2 /im is shown in figure 2 for each 
ion type. A simplified geometry is assumed in which 
the channel length and width and the collection 
length (fig. 3) are taken as equal to the feature size. 
The E, Z plane can be divided into regions for which 

A Q(E) > Q c (8) 

The value of A Q(E) depends on the feature size L. 
(See eq. (1).) The ion mass for each value of Z was 
taken as the natural mass in arriving at the contour 
of constant A Q shown in figure 4. The average 
recoil energies from the fragmentation of 16 O and 
28 Si produced by collision with high-energy protons 
(ref. 14) are also shown in figure 4. The importance 
of a given fragment type for a given feature size for 
the device may be judged from the average recoil 
energies from the fragmentation of ie O and 28 Si. 
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It is doubtful that any of the fragments produce 
upsets in the 4-/rm and larger devices (note that 
simplified geometries have been used). Also, the 
lighter fragments of Li, He, and H are not suspected 
for SEU events in this simple geometry and figure 4 
is applicable to incident cosmic-ray ions. 


Nuclear Recoil Transport 

The transport of the recoil fragments is described 
as follows: 


n-v-JL Sz ( E ) 


Mx,tt,E) = UE) 


( 9 ) 


Nuclear-Fragmentation Cross Section 

Although nuclear fragmentation has been under 
study for nearly 50 years, the absolute cross sections 
still stir some controversy. The experimental prob- 
lem was that the main-reaction products could be 
directly observed only in recent years and even now 
only in rather sophisticated experiments. Rudstam 
(ref. 16) studied the systematics of nuclear fragmen- 
tation and supposed the fragment isotopes to be in 
a bell-shaped distribution about the nuclear stabil- 
ity line. Silberberg and Tsao (ref. 12) continued the 
Rudstam parametric approach and added many cor- 
rection factors as new experimental evidence became 
available. 

Concurrently, Monte Carlo simulation of the Ser- 
ber model (ref. 17) and final decay through com- 
pound nuclear models showed some success (refs. 7 
and 18). Even so, Monte Carlo simulation shows 
little success in predicting fragments whose mass is 
small compared with the original target nuclear mass 
(refs. 13 and 19). Of the various models for nucleon- 
induced fragmentation in 28 Si, the model of Silber- 
berg and Tsao (ref. 12) is probably the most reliable. 
The main limitation of their model is that only inclu- 
sive cross sections are predicted; particle correlations 
could prove important in predicting SEU. 

Measurements of 27 A1 fragmentation in pro- 
ton beams have been made by Kwiatkowski et al. 
(ref. 19). These experiments are compared in figure 5 
with the Monte Carlo results of OMNI and MECC-7. 
Also shown are the results from reference 16; gener- 
ally, these results appear to be within a factor of 2 
of the experiment. The model in reference 12 is the 
only model which predicts significant contributions 
in the important range below A — 12. 

The spectrum of average recoil energy is calcu- 
lated using the formalism of reference 14 and the 
Silberberg-Tsao cross sections and is shown for com- 
parison with the spectrum according to the Bertini 
cross sections in figure 6. The Bertini cross sec- 
tions are greatly underestimated above 8 MeV and 
greatly overestimated below 3 MeV. The Bertini re- 
sults are typical for currently available intranuclear- 
cascade models. Experimental evidence indicates 
that even the Silberberg-Tsao values are too small 
above 6 MeV (ref. 19). 


where (f> z (x, S7, E) is the ion flux at x moving in 
direction Q, with energy E and where Cz(E) is the ion- 
source density assumed to be isotropic and uniformly 
distributed through the media. The solution to 
equation (9) is in a closed region bounded by the 
surface T subject to the boundary condition 

<j> z {r, f2, E) = F z (fl, E) (n • S7 < 0) (10) 

where n is the outward-directed normal of the sur- 
face T. The solution is found by using the method of 
characteristics (refs. 20 and 21) as 

0 2 (x,J2,£) = |^0 2 (r,ft,E 6 ) 

+ S Z (E) J E (11) 

where T is the point on the boundary determined by 
projecting x along in direction O and 

E b = R z 1 [R z (E) + b] (12) 

where 

b = n-(x-r) (13) 

Equation (11) may be used to evaluate the spectrum 
of particles leaving the region that can be related 
to the spectrum of energy deposited in the media. 
An isolated sheet of silicon of thickness a, which 
is obviously similar to the McNulty surface-barrier 
detectors, is considered. The inward-directed flux 
at the boundary is then zero. We first consider a 
monoenergetic ion source 

Cz(E) = ^ S(E - E ') (14) 

for which 

a & f 1 (E E ' "C Eh) 

Mx,n,E)= ~ b) 

47 tS z (E) ( o (Otherwise) 

( 15 ) 

where <r z is the silicon-fragmentation cross section 
and (f) is the flux of initiating energetic particles. 
The spectrum of ions leaving the sheet (ignoring edge 
effects) is 


3 



( 16 ) 


dfz 

dE 


= A-kA 



E) dfj 


A<T z <f> 
2 S Z (E) 


[R z {E')-R z {E)f 

1 

0 


(0 <E<R~ l [R Z {E' -a]) 

{R- z l [R z {E')-a]<E<E') 
(E f < E) 


where A is the area of the sheet and n is the cosine of the colatitude with respect to the local surface normal. 
The total number of escaping particles is found by integrating the spectrum given by equation (16) and is 


N e = Aaa z (f> 


i a 

1 2R Z {E')\ 

Rz(E') 

2a 


(a < R Z (E ' )) ) 

(a > R Z (E")) j 


From equation (17), the total number of ions which stop in the sheet is 


N s - Aaa z (f> < 


[2 Rz(E')} 


(a < R Z (E')) ' 
(a > R Z (E')) / 


(17) 


(18) 


Obviously, an ion produced with energy E' which leaves the sheet with energy E suffered an energy loss e to 
the sheet given by 


e = E f — E 


(19) 


which we use to find the energy-loss spectrum as 


dfz6 

de 


dE 


+ N S 6(E' - e) 

E—E'—e 


(20) 


Considering that equation (20) is the energy deposition in a sheet of area A and thickness a as the result of a 
monoenergetic volumetric source, the response to any arbitrary spectral source can be found by superposition. 


Fragmentation Energy-Loss Spectra 


The fragmentation-source energy distribution (normalized to unity) is given as 


p(£,) = /5 exp teS) 


where 3Eq is the mean-fragment energy and is given by Wilson et al. (ref. 14) based on previous work by 
Goldhaber (ref. 22). 
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The energy-loss spectrum is found by using equations (20) and (21) as 


dF = 7°° df z6 
de J € de 

= f 30 dfzs 
Jo de 


p(E') dE' 

E=E'-e 

p(E + e) dE 

E=E'+e 


The contribution from stopping ions is readily evaluated to give 


( 22 ) 


dF 

de 


Ns{c)p(e) + 



dfz6 

de 


p(E -(- e) dE 

E'=E+e 


(23) 


where the second term of equation (23) requires more attention. 

The energy-degradation function in the integral of equation (23) is given by equation (16). It is not clear 
how the integral in equation (23) is to be evaluated. As an approximate evaluation, the energy-degradation 
function is approximated by two or three line segments as shown in equations (24) to (27). 

If Rz(e) > 2a, then 


df z 

dE 


E'=E+e 


Acr z (f) 
2 S Z (E) 


RUe) 


1 _ a 2 

? J8JT) 


Rz(E) 

R^E~2 1 


where E 2 is the solution of 
and E\ is the solution of 

In the event that Rz(c) < 2a, then 


dfz 

Ao z 4> 

f a2 4 . 

1 “ 2 
Wtn 

dE 

E'=E+£ - ] 

1 

1 


1 3[R Z (E)-R Z (E 2 )} 

4 + A[R z {Ei)-R z {E 2 )\ 

1 

Rz( E 2 ) = Rz{E 2 + e) - 2 a 
Rz(Ei) — Rz(E 1 + e) — a 


Rz(E) 

RjfEi T 


(0 < E < E 2 ) 

(E 2 < E < E\) 

{E\ < E < 00 ) , 


(0 < E < El) 

> 

(-Ei < E < 00 ) 


(24) 


(25) 

(26) 


(27) 


with the understanding that E\ is zero if Rz(t) 

< a. 

The second term of equation (23) 

is divided into three 

subintervals as follows: 





h(<) = 

fE 2 

dh 

de 

p(E + e) dE 

(28) 


Jo 

E'=E+e 

h(t) = 

rEi 

dfz 

de 

p(E -1- e) dE 

(29) 


Je 2 

E'=E+c 

73(f) = 

roc 

dfz 

de 

p(E + e) dE 

(30) 


Jei 

E'=E+e 


First, I\{e) is zero unless Rz(c) > 2a, for which 
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h(e) = 


Aa z (f) J a‘ 


#zO) 


P(E 2 ,e) + 


R z( e ) 


Q(E 2,<0 
R z{E 2) 


/2(e) = 


Ac z (f> 


(j 


R z{E 2) 


4 [i? z (^i) - ^(£2)] 


[P(£7i,c)-P(£? 2 ,e)] 


+ 


Aa z <j> 3 Q(E\,e) - Q{E 2: e) 
2 A[R z {E\) - Rz{E 2 )\ 


A(T z <t> f°° p(E + e) 




If a < Rz{t) < 2a, then E 2 and /i(e) are zero and 

Aa z (j> ( a 2 


h(e) = 


R%M) 


P(E i,e) + 


1 - 


R 2 M 


Q(£i,<) '| 

Rz(Ei) J 


When Rz{t) < a, then E\ = E 2 = 0, so that /i(e) and I 2 (e) both vanish and 

Acr z (j) f°° p{E + e) 


^ = -^rk r -fzw dE 

In equations (31), (32), and (34), P and Q are given by 

* Ei p(E + e) 


P(Ei,e)= f 

Jo 


0 S Z (E) 


dE 


Q(E h e)= [ 

Jo 


Ei R z (E)p(E + e) 
S Z (E ) 


dE 


The integral of equation (36) may be approximated for values of Ei < by 

d/it ^ _ Rz( E o) N fl Ei \ 

P(Ei,e)K—j=-p(e) 1 

where 7 is an incomplete gamma function (ref. 23). For larger values of Ei (j^c < E{ < 4e^, the 
be taken as 

R z(Eo)_ nU ,f IJ1 + 1 n 

1 r» / \ r% ’ 


p {Et ,e ) « {i 7 (1, 4 ) + 17 ( 2 . 


+ 




Whenever Ei > 4e, the integral is approximately 

{r (i sk ) + r (I- i) 




(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 

integral may 

(39) 

(40) 
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The integral in equation (37) may be approximated 

by 


Q(Ei,e) « R \^ o) { C(E t + e) - C(e)] (41) 

LhjQ 

where C(e) is the integral spectrum as follows: 

C(E) = [ E p(E') dE' (42) 

Jo 

A useful check on the approximations involved is the 
strict requirement 


Aa z (f) f°° p(E + e) 




The total absorption spectrum is then 
dF 


de 


— Ng{e)p{e) + Ji(e) + h{t) + hi e ) (44) 


and is shown in figure 7 for detector thicknesses of 
1 to 5 pm with Eq — 3.5 MeV. Similar results are 
shown in figure 8 for detector thicknesses of 50 to 
200 pm. In comparing figures 8 and 9, it is shown 
that the energy-loss spectrum is approaching the 
fragment-production spectrum as a becomes larger. 
The normalization is always 


r°° dF_ 

Jo de 


= 1 


(45) 


which is satisfied by numerical evaluation to within 
2 percent. 


Results 

Typical fragmentation cross sections calculated 
using the Silberberg-Tsao model are shown in ta- 
ble 1 for 125-MeV protons. The values of Eq are 
taken from reference 14. The calculated response 
of the 2.5-jum detector is shown in figure 10; these 
values should be compared with the experiments of 
McNulty and the values according to the Monte 
Carlo code of the McNulty group, which are also 
shown in figure 10. The peak value at zero energy 
is fixed by the total reaction cross section and total 
proton flux. It appears that the total reaction cross 
section of the McNulty code is too small. Otherwise, 
the present theory and the Monte Carlo code show 
nearly equivalent agreement with the experiments. 
Similar comments apply to the 4.2-/im detector re- 
sponse (fig. 11) with one exception. The energetic 
events above 20 MeV observed in experiments are 
well represented by the present theory but not by the 


Monte Carlo code as expected by observing figure 6. 
This high-energy agreement between theory and ex- 
periment is observed for the 24.1-/im detector, but 
the Monte Carlo code again fails to predict the high- 
energy events as shown in figure 12. The improved 
model of the present work is again clearly displayed 
for the 158-MeV experiments of McNulty et al. as 
shown in figures 13 and 14. 

The inability of the Monte Carlo code to predict 
the most energetic fragments could be a serious lim- 
itation in predicting SEU in some devices. Although 
the Silberberg-Tsao cross sections for proton-induced 
reactions are not in complete agreement with some 
recent cross-section measurements, they still provide 
improved ability over Monte Carlo models. The 
methods of analysis used in this paper will be ap- 
plied to specific-device geometries in the near future. 

Concluding Remarks 

Although the intranuclear-cascade code provides 
a very detailed model of the nuclear reactions of high- 
energy protons, the inaccuracies of current cascade 
and compound nuclear-decay models limit the use- 
fulness of such models in single-event upset stud- 
ies. The improved predictions of the current semi- 
empirical cross sections of the Silberberg-Tsao model 
seem adequate for microelectronic studies. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
September 7, 1989 
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Table 1. Cross-Section Parameters for Fragmentation of 
28 Si by 125-MeV Protons 









Charge collected, pC 



L, jim 

Figure 1. Critical charge as a function of feature size in several device types. 
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Ion energy, MeV 

Figure 2. Charge collected as a function of ion energy with a collection length of 2 ^m. 
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Figure 3. Cross section of bulk CMOS technology. 



Figure 4. Critical energy as a function of ion charge for several feature sizes. Average recoil energies of 
fragments of silicon and oxygen are superimposed. 
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with experimental measurements. 
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Figure 7. Total absorption spectrum for surface- barrier detector of 1 to 5 fim and Eq = 3.5 MeV. 



Figure 8. Total absorption spectrum for surface- barrier detector of 50 to 200 // m and Eq = 3.5 MeV. 
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Figure 9. Energy-loss spectrum for Eq — 3.5 MeV. 
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Integral number 



Figure 11. Response of a 4.2-/im surface-barrier detector to 125-MeV protons (2.14 x 10 8 protons). 
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Figure 13. Response of a 2.5-/rm 


surface-barrier detector to 158-MeV protons (3.9 x 10® protons) 



Integral number 



Figure 14. Response of an 8.7-/rm surface-barrier detector to 158-MeV protons (3.9 x 10 8 protons). 
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